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Subsurface grid element 
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Modeling of basin of fractal 
geometry using Voronoi 
tessellation 



Input to the program 

1. Boundary limits of the region 

2. Voronoi centers with in the boundary 

3. The value of p for L p norm 

4. Density values corresponding to each region (could be same also) 

5. Number of grid nodes in x and y direction 

6. Grid node spacing in x and y direction. 

7. Depth of the sub-surface region 



Generation of fractal sub-surface and density 
assignment to each region. 
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Laid grid of the specified specifications at the 
surface over the tessellated region of interest. 



Computation of gravity response due to each 
subsurface polygonal area of different physical 
property in tessellated region at each nodes of the 
grid at the surface. 



Cumulative gravity response at each nodes of the 
grid laid at the surface. 
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Increase the depth value and repeat all steps 
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Cumulative sum of the gravity anomaly with 
respect to depth at each nodes of the grid at the 
surface. 
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